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Abstract 

The spatio-temporal dynamics of three interacting species, two preys and one predator, in the 
presence of two different kinds of noise sources is studied. To describe the spatial distributions of 
the species we use a model based on Lotka-Volterra equations. A correlated dichotomous noise 
acts on (3, the interaction parameter between the two preys, and a multiplicative white noise 
affects directly the dynamics of each one of the three species. We study the time behaviour of 
the three species in single site for different values of the multiplicative noise intensity, finding 
noise-induced oscillations of the three species densities with an anticorrelated behaviour of the two 
preys. Afterwards, by considering a spatially extended system formed by a two-dimensional lattice 
with N sites and applying a mean field approach, we get the corresponding moment equations 
in Gaussian approximation. Within this formalism we obtain the time behaviour of the first and 
second order moments for different values of multiplicative noise intensity, with I3{t) subject to the 
same dichotomous noise source. Finally, we compare our results with those obtained by using a 
coupled map lattice model, consisting of a time discrete version of the Lotka-Volterra equations. 
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I. INTRODUCTION 

Noise is not generally detrimental to biological systems but can be employed to generate 
genotypic, plienotypic, and behavioral diversity 3, 4|. Real ecosystems are affected by 
the presence of noise sources which consist of random variability of environmental parame- 
ters, such as temperature, food availability, general conditions which can favour or thwart 
the increase of some biological species. This randomly fluctuating behaviour can be mod- 
eled by Gaussian noise sources, which influence, through a multiplicative interaction, the 
system dynamics. Multiplicative noise often causes the appearance of fluctuating barriers 
or processes of anomalous diffusion and has been investigated in the context of population 
growth and extinction Q, [s, ^, Q, y, [o], Q, 11, 12, 13, Q, Q- In this paper we study the 



time evolution of three interacting species, two preys, x and y, and one predator, z. The 
interaction between the two preys is symmetric and it is given by the parameter /3. We study 
the ecosystem dynamics, described by generalized Lotka-Volterra equations, in the presence 
of two different kinds of noise sources: (i) a dichotomous noise acting on the (3 parameter, 
(ii) three external sources, modeled as independent multiplicative Gaussian noises, which 
act directly on the three species. First we consider the deterministic dynamics of the system 
in a single site and we get the time behaviour of x, y and z, by analyzing the stability of the 
ecosystem with different constant values of the interaction parameter /3, which correspond to 
a coexistence regime {Pdown < 1) or to an exclusion regime {Pup > !)• Then we consider the 
interaction parameter f3 varying dichotomously between these two values. In this condition 
we study the time behaviour of the species concentrations x, y and z for different levels of the 
multiplicative noise intensity. We find noise-induced oscillations and strong anticorrelations 
between the preys. Afterwards we take into account the spatial version of our ecosystem, 
considering a two-dimensional domain formed by sites and adding a diffusion term in 
the L-V equations. By using a mean field approach, we obtain the corresponding moment 
equations in Gaussian approximation. We find that, for j3 varying dichotomously, the 1** or- 
der moments of the three species concentrations are independent on the multiplicative noise 
intensity. On the other hand, the behavior of the 2"'^ order moments is strongly affected 
by the presence of external noise sources. In particular we find that the time behavior is 
anticorrelated for the species densities of the two preys, and correlated between the predator 
and the total density of the two preys. Finally we get the time behavior of the 1** and 2"'^ 
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order moments using a coupled map lattice (CML) model [l6j and we compare these results 
with those previously obtained within the mean field approach. In view of an application on 
real systems, the results obtained could be useful to explain experimental data, reproducing 

nnnn 

the behaviour of natural ecosystems [5|, la, [171, 118[ . 



II. THE MODEL 

Our system is described by a time evolution model of Lotka-Volterra equations, within 
the Ito scheme, with diffusive terms in a spatial lattice consisting of sites 

Xij = X Xij (1 - z/ Xij - (3 Hi J - a Zij) + Xij Vo^^ij + D{< x > -Xij) (1) 
yi,j = A Vij {l-u yij - 13 Xij - a Zij) + y,j ^fo'yUj + D{< y > -yij) (2) 
Zi,j = Xz Zij [-1 + 7 {xij + yij)] + Zij -Jolilj + -D(< z > -Zij), (3) 

where the dot indicates the time derivative. The variables Xij, yij and functions 
of the time t, and denote the densities, respectively, of the two preys and the predator in 
the lattice site («, j). A and A^ are scale factors, z/ is the growth rate for the two preys, D 
is the diffusion coefficient, and <x>,<y>,<z> indicate the spatial mean, performed 
on the whole lattice, of the three species densities. The coefficient (3 is the interaction 
parameter between the two preys. The coefficients a and 7 account for the interaction 
between preys and predator, (.fjit), C,ij(t), Hjit) are statistically independent Gaussian 
white noises with zero mean and unit variance, and they model the interaction between 
species and environment. Finally, a^, Cy, are the intensities of the three sources of 
Gaussian white noise. 



A. Single site dynamics 



1. Stability analysis and dynamical regimes 

Depending on the value of the interaction parameter, coexistence or exclusion regimes 
take place. In the absence both of multiplicative noise {a^ = = o"^ = 0) and diffusion 
terms {D = 0), Eqs. ([I])-© describe the deterministic dynamics of a single site ecosystem. 
In these conditions, for the generic site of lattice the stationary values of the three species 
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densities are given by 



^stat ^ yStat ^ }_ 



= h^SllA, (5) 

where the indices i,j where suppressed. From Eqs. ([I])-© one can see that the two prey 
densities have stationary values that are independent on the interaction parameter (3. Con- 
versely, the stationary value of the predator density is connected with the value of (3. This 
indicates that the interaction parameter between the two preys determines the coexistence 
or exclusion regimes for the whole system, affecting the stationary value 2;'**"*. From Eq. ([5]) 
the survivance condition for the predator is 2**'^* > 0, which allows to get the coexistence 
condition for the three species as a function of (3 

/3 < 27 - z/. (6) 

The inequality ([6]) indicates that the system is characterized by two stationary states, which 
become stable or unstable depending on the values that (3, 7 and v take on. In particular, 
when the condition ([6]) is satisfied, the stable state is represented by the coexistence of 
the three species. Otherwise, after a transient, the predator tends to disappear (inequality 
([6]) doesn't hold anymore) and we get a system formed by two competing species, whose 
coexistence/exclusion conditions depend directly on the value of the parameter (3 jl?!. [l9 
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22I ]. In this sense, the predator plays a regulatory role for the dynamics of the two preys, 
whose reciprocal behavior is mediated by the interaction parameter (3 through the presence 
of the species z. We calculate the numerical solutions for single site dynamics, setting in 
Eqs. (HD - © A = 3, = 0.06, z/ = 1, a = 0.02, 7=1, with two different values l3down = 0.94 
and Pup = 1-04 of the interaction parameter /? and initial conditions x(0) = y{0) = 0.1, z{Q) 
= 2.0. The values of multiplicative noise intensity are the same for the three species, that 
is cr = = o-y = (Jz- In Fig. [T]we show the time series of the three species in coexistence 
iP = Pdown) and exclusion (/5 = Pup) regimes, for a = and cr = 10^"*^^. We note that, 
when the system is subject to deterministic dynamics, the coexistence regime causes, after 
a transient, the three species to reach the equilibrium values, x'^*"* = y'**"* = 0.5, 2;**"* = 1.5, 
obtained from Eqs. (H]) - (El) using u = 1, a = 0.02, 7 = 1,/? = Pdown = 0.94 (see Fig. [IK). In 
deterministic exclusion regime the predator tends very slowly to vanish. However, the two 
prey densities reach the stationary values, remaining constant (Fig. Wp)- In this case, the 
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FIG. 1: Time evolution of the three species. Deterministic dynamics in (a) coexistence and (b) 
exclusion regime. Stochastic dynamics, for a = 10^^^ in (c) coexistence and (d) exclusion regime. 
Values of the parameters and initial conditions are A = 3, = 0.06, u = 1, a = 0.02, 7 = 1, x(0) 
= y{0) = 0.1, z{0) = 2.0. 

stationary values correspond to an unstable equilibrium point. In fact, in the presence of 
a small level of multiplicative noise, the symmetry, due to the parameter values and initial 
conditions used in our simulations, is broken and one of the two preys prevails, displacing the 
other one (Fig. [T]d), according to the previously obtained results 2l|. Finally we note that no 
significative modifications occur, with respect to the deterministic case, when a small level 
of noise is present in coexistence regime (see Fig. [It). This obviously depends on the fact 
that, for P = Pdown, the system occupies a stable equilibrium point, which is maintained also 
in the presence of low levels of multiplicative noise. However, environmental perturbations, 
due to the presence both of deterministic and random fluctuations of biological and physical 
variables, such as the temperature, affect the dynamics of the species. These external forces 
can modify the behaviour of the populations, either introducing multiplicative noise sources 
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FIG. 2: Time evolution of the interaction parameter (3{t) with initial value /3(0) = 1.04 and delay 
time Td = 435. The interaction parameter I3{t) switches quasi-periodically between Pdown = 0.94 
and f3up = 1.04. The values of the other parameters are A = 9.0, uj/{2'k) = 10~^, xo = 2 • 10^^. 

which act directly on the species or affecting the dynamics of the interaction parameter 
f3. In fact, the environmental variations can cause the system dynamics to change between 
coexistence (/5 < 27 — u) and exclusion (/? > 27 — u) regimes. This dynamical behavior 
can be described by considering that the interaction parameter j3{t) is a stochastic process 
driven by a dichotomous noise, whose jump rate is given by 

X{t) = { (7) 

Xo(l +^|cOSCjt|), At > Trf . 

where At is the time interval between two consecutive switches, and Td is the delay be- 
tween two jumps, that is the time interval after a switch, before another jump can occur. In 
Eq. ([7]), A and u = {2n)/T are respectively amplitude and angular frequency of the periodic 
term, and Xo is the jump rate in the absence of periodic term. This causes P(t) to jump 
between two values, Pdown < 27 — z/ and Pup > 27 — z/. According to the condition ([2D, 
these values determine the two possible dynamical regimes (coexistence or exclusion) of the 
deterministic Lotka-Volterra's model for three interacting species. For given values of the 
parameters A, uj and xo the switching time between the two levels of /9(t) depends on r^. 



Applying a procedure analogous to that followed for the two-species case [19|, we set A = 9.0, 
u/ {2tt) = 10~^, xo = 2-10~^, obtaining the time series of Pit) for = 435, with Pdown = 0.94 
and Pup = 1.04. The results, shown in Fig. [21 indicate the presence of a synchronization 
between the jumps and the periodicity of the rate x{t)- Fo^' ^ system formed by two com- 
peting species this causes a quasi-periodical time behavior of the two populations, which 



can be considered as a signature of the stochastic resonance phenomenon 2J] in population 



dynamics 
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|22|. Therefore we fix the delay at the value td = 435, which determines an 
oscillating dynamical regime. In these conditions, f3{t) switches quasi-periodically between 
Pdown and Pup (see Fig. [2]), causing the system to be alternatively subject to the coexistence 



and exclusion regimes. 



2. Time behaviour of the species in a single site 

In this section we analyze the time behaviour of three interacting species in a single site 
of the lattice. From Eqs. ([I])-®, by setting D = we get 

X = \x {1 — ux — Py — az) + X y/ox^^(t) (8) 
y = Xyil-uy-(3x-az)+y^e{t) (9) 
z = X,z{-l+^x + jy) + z,/r,C{t), (10) 

where the indices i,j where suppressed. 

By choosing /3(0) = 1.04 and = 435, we obtain for i3{t) the time behaviour shown 
in Fig. [21 We analyze the time evolution of the species densities by numerical simulation 
of Eqs. (l8l)- (fT0l) . The time series of x, y and z are obtained for different values of the 
multiplicative noise intensity, namely a = 0, 10~^^, 10~^, 10"'^. The values of the other 
parameters are the same used in the previous section, that is A = 3, = 0.06, z/ = 1, 
a = 0.02, 7 = 1, Pdown = 0.94, Pup = 1.04. The initial values of the species densities are x{0) 
= y{0) = 0.1, z{0) = 2.0. In Fig. [31 where the results are reported, the time series of x(t), 
y{t) (preys) and z(t) (predator) show correlated behaviour in the absence of noise (panel 
a). In the presence of noise intensity an anticorrelated oscillating behaviour of x{t) and 
y{t) appears (see panels (b)-(d)). Moreover we note that, for all the values of multiplicative 
noise intensity, the two prey densities oscillate, with the frequency of the external driving 
force, around the stationary values, x''*"* = y''*"* = 0.5. We observe that the predator 
density show an oscillating behaviour, with the same frequency, around a value smaller than 
^stat _ g_ However, the oscillations of z{t) are characterized by a larger amplitude with 
respect to < x{t) > and < y{t) >. This behaviour is connected with the different effect 
that the alternating regime (exclusion/coexistence) produces on preys and predator. In 
fact, the quasi-periodical behaviour of P{t) affects directly the dynamics of the predator (see 
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FIG. 3: Time evolution of the three species densities in a single site of the lattice. The values of 
the multiplicative noise intensity are: (a) a = 0, (b) a = 10~^^, (c) a = 10^^, (d) a = 10^'^. Here 
A = 3, A2 = 0.06, u = 1, a = 0.02, 7 = 1. The values of the other parameters are the same of 
Fig. [2j The initial values of the species densities are x(0) = y(0) = 0.1, 2;(0) = 2.0. The time 
series of x{t), y{t) (preys) and z{t) (predator) show a correlated behaviour in the absence of noise 
(panel a). In the presence of the noise (panels (b)-(d)) an anticorrelated behaviour of x{t) and y{t) 
appears. 

Eq. [5]), causing a decrease of the mean value of z during the exclusion regime. Conversely, 
in coexistence regime the two preys maintain a constant value (see Eq. (jlj)) going towards 
an anticorrelated regime for (3{t) = (3up- In this last condition the two preys are subject to a 
pure competitive dynamics, recovering the behaviour observed in a system of two competing 



species 



2l| 
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B. Spatially extended system: Mean field approach 

In this section we analyze the time behaviour of three interacting species in a spatially 
extended system by using a mean field approach. The system dynamics is described by 
Eqs. ([I])-© in the presence of the diffusive term {D ^ 0). In order to use a mean field 
approach we derive the moment equations for this system. Assuming — > oo, we write 
Eqs. ([I])-® in a mean field form 

X = fxix, y, z) + v^^,;(a;) Cif) + D{< x> -x), (11) 
y = fy{x,y,z) + ^jgy{y)^yit)+D{<y> -y), (12) 
z = Mx, y, z) + ^ygM C{t) + D{<z> -z), (13) 

where < x >, < y > and < z > are average values on the spatial lattice considered (ensemble 
averages in the thermodynamic limit) and we set fx{x, y, z) = \x{l—vx—(3y—az)^ Qxi^) = x, 
fy{x, y, z) = Xy{l - uy - /3x- az), gy{y) = y, f^{x, y, z) = Xzz[-1 + 7(x + y)], g^{z) = z. 
By site averaging Eqs. (fTT!) - (fT3l) . we obtain 

< a; > = <f^{x,y,z)>, <y> = <fy{x,y,z)>, <z> = <f^{x,y,z)>. 

(14) 

By expanding the functions fx{x,y,z), gx{x), fy{x,y,z), gy{y), fz{x,y,z), g^iz) around the 
1^* order moments < x{t) >, < y{t) > and < z{t) >, we get an infinite set of simultaneous 
ordinary differential equations for all the moments 25j. To truncate this set we apply a 
Gaussian approximation, for which the cumulants above the 2"'^ order vanish. Therefore we 
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obtain 

< X > = X<x>(l — i'<x> —[3 < y > —a < z >) — A(z/ /i2oo + P A'-iio + A'-ioi) (15) 

<y> = \<y> {l-v <y> -(3 < x > -a < z >) - \{v /io2o + (3 A^iio + a mi) (16) 

<z> = A^<2;>(-l + 7<x> +7 <?/>) + A^7(/iioi + /ion) (17) 

/i2oo = 2A(1 -2iy < X > -(3 < y > -a < z >)/i2oo 

- 2A < X > iPfiuo + a/iioi) + 2a^{fi2oo+ < x >^) - 2L)/i2oo (18) 
/io2o = 2A(1 -2z/ <y>-f3 < x > -a < z >)/io2o 

- 2A < y > (/5/iiio + a/ioii) + 2(Ty(/io2o+ < y >^) - 2L)/io2o (19) 
/ioo2 = 2A^(-1 + 7 <a;>+7 <?/ >)/xoo2 

+ 2A^7 <z> (/iioi + /ion) + 2(T^(/ioo2+ < z >^) - 2D/ioo2 (20) 

/iiio = X[2 - 2iy {< X > + < y >) - /3 {< X > + < y >) - 2a < z >]/fiio 

- A/5(< X > /io2o+ < y > /i2oo) - Aa(< a; > /ion+ < y > /^loi) - 2D/iiio (21) 
/iioi = A(l - 2z/ < X > -/5 <y> -a < z >)/iioi + A^(-l + 7 <x>+7 <y >)/iioi 

- A < X > (a /ioo2 + /5 /ion) + A^7 < 2; > {fiuo + /i2oo) - 2L'/iioi (22) 
/ion = A(l - 2z/ < y > -/3 < X > -a < z >)/ion + + 7 <a;>+7 <?/ >)/ion 

- A < ?/ > (a /ioo2 + P /iioi) + A^7 < 2; > (/ino + /io2o) - 2L'/ion- (23) 

where /t2oo, /i020, /ioo2, /^iio, /iioi, /ion are the 2"'^ order central moments defined on the 
lattice 

/^200(i) = < X^ > - < X >^ /i02o(^) = < i/^ > - < 2/ >^ /i002(^) = < 2^ > - < ^ >^ 

/ino(^) = < xy > — < X >< y >, /iioi(^) = < xz > — < x >< z >, /ion(^) = < yz > — < y >< z 

In order to get the dynamics of the three species we analyze the time evolution of the 1^* 
and 2"''^ order moments according to Eqs. (fT5l) - (!23l) . As initial conditions we consider each 
species uniformly distributed on the spatial domain, that is we set < x(0) > = < y{0) > 
= 0.1, < z{0) > = 2.0, /t2oo(0) = /io2o(0) = /ioo2(0) = /ino(O) = /iioi(O) = /ion(O) = 0. 
Therefore, from Eqs. ffT^ - fl25]) we get, in the deterministic case, the stationary values for 
< X >, < y > and < z > 

<a;>^*«* = <y = — , < z = ~ ^ ''^ (25) 

27 2q;7 
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Using for the parameters the same values of the single site analysis, we obtain < 
^ -ystat _ ^ y -ystat _ ^ ^ ^stat _ ^^^so fix the delay time at 

the same value = 435 used in the single site case. Finally, by numerical integration of 
Eqs. ffT^ - fl2^ . setting D = 10^^, we get the time series of the I''* and 2"°' order moments for 
the following values of multiplicative noise intensity a = 0, 10^^^, 10^^, 10^^. The results are 
reported in Figs. HJ [51 Here we note that, after a transient, the mean values of the two prey 
densities (see panels a and d of Figs. HI [5]) oscillate around the stationary values. The oscil- 
lations are connected with the presence of two stable equilibrium points. For /? = Pdown < 1 
the stable equilibrium is given by the contemporary presence of the three species (coexistence 
regime). Conversely, for (3 = f3up > 1 the system goes towards a new equilibrium point, with 
the predator tending to disappear (exclusion regime). In the presence of a dynamical regime 
(the system switches periodically from coexistence to exclusion), we observe the appearance 
of correlated oscillations in the time series of < x(t) >, < y(t) > and < z{t) >. In particular, 
we note that < z(t) > is subject to oscillations occurring around a value smaller than the 
stationary one (< z >'^*'^*= 1.5) and characterized by a larger amplitude with respect to 
< x{t) > and < y(t) >. This behaviour is analogous to that observed in the case of single 
site dynamics. In the absence of noise (top of Fig. H]), the time series of < x{t) >, < y{t) > 
and < z(t) > (panel a), fi2oo(t), /io2o(^), /ioo2(t) (panel b) and /iiio(t), /iioi(^), /ioii(i) (panel 
c) are completely overlapped and each species maintains a homogeneous distribution over 
the lattice, that is all the 2"'' order moments remain equal to zero. For cr = 10^^^ (bottom 
of Fig. H]) no changes are observed in the behaviour of the mean values (see panel d), and 
the variances of the three species show correlated oscillations (panel e). In panel f, fiuo 
oscillates taking on only negative values. This indicates that the spatial distributions in the 
lattice will be characterized by the presence of regions where species x or species y prevails. 
The two preys will be distributed therefore in non-overlapping spatial patterns. This picture 
is in agreement with previous results obtained with a different model 231]. Conversely, /xioi 
and yUoii are always zero (see panel f of Fig. Hj). This behaviour indicates that the predator 
is uncorrelated with the density of each prey: the species z tends to occupy indifferently 
the sites where x oi y prevails (see the time behaviour of /ioo2 in panel e of Fig. H]), but is 
correlated with the total prey density (a global increase of food availability improves the 
life conditions of the predator) . This explains why the variance of the predator shows small 
oscillations. On the other hand, when exclusion regime takes place, the two preys tend to 
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FIG. 4: Time evolution of the 1'^* and 2 order moments in the mean field approach. The values 
of the multiplicative noise intensity are: a = 0, 10^^^ from top to bottom. In the absence of noise 
the time series of < x{t) >, < y{t) > (panel a), fi2oo, A^020, Ai002 (panel b) and /xno, ^loi, /^oii 
(panel c) are completely overlapped. The predator (mean value of species z) shows a behaviour 
correlated with those of both preys (mean values of species x and y). For a = 10~^^, no changes 
are observed in the behaviour of the mean values (panel d), the variances of the two preys oscillate 
overlapping each other and a correlation is observed with the variance of the species z (panel e), 
the covariance of the two preys, //no, oscillates taking on only negative values (the two preys 
are anticorrelated each other), while /iioi and /ion are always zero (panel f). The initial values 
of the moments are < x(0) > = < y(0) > = 0.1, < 2;(0) > = 2.0, /i20o(0)=/io2o(0)=/ioo2(0) = 
/^iio(0)=/^iio(0) = 0=/ioii(0)=0. The diffusion coefficient is D = W~^. The values of the other 
parameters are the same used in Fig. [3l 

occupy different sites, "spreading out" in the spatial domain and causing an increase of their 
variances (see panel e of Fig. Hj) with a stronger anticorrelation (see the behaviour of /ino 
in panel f of Fig. Hj). Finally we note that the amplitude of the oscillations both of the 
variances and covariances increases as a function of the noise intensity: in particular they 
have the same order of magnitude of a (see panels b, c, e, f, in Figs. HJ [5]). In fact, for higher 
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FIG. 5: Time evolution of the 1'^* and 2"'^ order moments. The values of the multiplicative noise 
intensity are: 10~^, 10^'^ from top to bottom. No changes are observed in the time behaviour of 
<j;> = <y> = <2;> (see panels a and d) for both values of the noise intensity. An increase 
in the amplitude of oscillations, as a function of the noise intensity, appears both in the variances 
of the two predators, fi200, Ai020) (see panels b and e) and in the covariance of the two preys, Huo, 
(see panels c and f). The values both of initial conditions and parameters are the same used in 
Fig. El 



levels of multiplicative noise {a = 10~^, 10^^) the amplitude of the oscillations increases 
and the periodical anticorrelated behaviour between the two preys becomes more evident. 
Conversely, no modifications appear in the time series of the mean values as a function of 
the multiplicative noise intensity (see panels a, d in Figs. HIED. 

Even if it is related to a very different mechanism, this behavior is similar to the stochas- 
tic resonance effect produced in population dynamics, when the interaction parameter is 
subjected to an oscillating bistable potential in the presence of additive noise 



21 



22 . We 



note that in the absence of external noise (cr = 0) both pop ulations coexist and the species 
densities oscillate in phase around their stationary value 211]. This occurs identically in each 
site of the spatial lattice (single site dynamics). The behavior of the mean values reproduces 
this situation. For cr 7^ 0, in the single site dynamics we observe anticorrelated oscillations of 
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X and y (preys). By site averaging these noise- induced oscillations (see Ref. [2l|]) we recover 
the average behavior obtained in the absence of noise. This spatial auto-averaging effect 
explains why the 1*** order moment behavior is independent on the external noise intensity, 
while the 2"^^ order moments give information on "spreading" and anticorrelation of the 
species densities in the spatial domain. 

III. COUPLED MAP LATTICE MODEL 

In this section we adopt a different approach to analyze the dynamics of the three species 
on the square lattice defined in Section [TTl We consider the time evolution of our system 



by using a coupled map lattice (CML) model IGj. In this formalism both correlated and 
anticorrelated spatial patterns of the three interacting species have been found j^]. Here 
we calculate the moments by using the CML model. By this approach, the dynamics of the 
spatial distributions of the three species is given by the following equations 

x^"^ = Xx^il - vx^ - - a 4^) + V^x^^f + D Y.i^f - x^\ (26) 

p 

= AygHl-^yi?-/3^"M?-«4?) + V^I/ffe^^ (27) 

p 

A.4?(-l + 7 + 7 + + D J2{z(;^ - z^), (28) 



Jn+i) 



where , ?/-"^ and z^j'' denote respectively the densities of prey x, prey y and predator 
z in the site {i,j) at the time step n. According to the notation used for the mean field 
approach. A, A^, u, (3, a, 7 and D represent the same quantities defined in Section [TTl 

Ci'f'' cire independent Gaussian white noise sources with zero mean and unit variance. 
The interaction parameter corresponds to the value of /5(t) taken at the time step n, 
according to Eq. ([7j). Here indicates the sum over the four nearest neighbours. 



A. Stationary states for the CML model 

Applying a procedure analogous to that used for Eqs. ([I])-(l3l), we consider Eqs. ( I26l) - (!28ll 
in the absence both of noise sources and diffusion terms (D = 0). In this conditions, for 
^ij-^^^ = yt^j^^^ = Vij^ -^ij^^"* ~ '^li'*' obtain the stationary values of the three species 
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densities for the generic site 



X 



stat 
CML 



stat 
^ CML 



1 

2^ 



^stat 



2^7 



(29) 



(30) 



where the indices were suppressed. As in the approach based on tlie use of differential 
equations, the stationary values of the two prey densities are independent on the interaction 
parameter /?, which is responsible for the two different dynamical regimes, coexistence or 
exclusion, and affects the dynamics of the whole system through its action on the stationary 
value 2;**"* . The existence condition for the predator 

CML -"^ 



^stat 



A3+1 



2a;7 



> 



(31) 



allows to get the following inequality for the interaction parameter (3 



V. 



(32) 



The inequality ( !32l) indicates that, according to the analysis performed in Section III Al 
the CML model is characterized by two stationary states that become stable or unstable 
depending on the values of the parameters. Comparing the inequalities (jSD and (15^ . we 
note that in the CML model the coexistence condition and the regulatory role, played by 
the predator on the dynamics of the two preys, depend also on the scale factors A and A^. 



B. Time series in the CML model 

In view of a comparison between mean field approach and CML model, we define the 
1^* and 2"°' order moments on the discrete lattice, at the time step n. The mean values, 

< X >("\ < y >("\ < z >("'\ given by 



<u> 



in) 



N 



{u = x,y,z) 



(33) 



represent the 1'^* order moments. The variances var^\ var^\ vari^^ defined as 



var. 



(n) 



N 



{u = x,y,z), 



(34) 
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and the covariances 

cov^2 = ^^^^ j;!^ in, w = x,y,z, u ^ w) (35) 

are the 2"'^ order central moments. In order to get, for the species densities, stationary 
values close to those obtained in the mean field approach (see Eqs. (l25l) ). we choose for all 
the parameters, except 7, the same values of Section [Tll Therefore, setting A = 3, = 0.06, 
z/ = 1, a = 0.02, 7 = 26.5, from Eqs. (!29l) . (!30|) . we calculate the stationary values for the 
densities of the two preys and predator in the coexistence regime (/3 = l3down = 0.94) 

< ^ >ctiL = <y >CML = 0.3 ; <z >%t^ = 1.0. (36) 

The CML model can be considered as a time discrete version of the Lotka-Volterra system, 
with time step At = 1. For the numerical integration of Eqs. ( |T5l) -(l23l) we used dt = 10~^, 
which is a suitable value to obtain convergence of the solution. Obviously, with these values 
of At and dt, the dynamics of the CML model results to be faster with respect to that 
obtained within the moment formalism. In particular, for (3 = (3up > 1, using the same 
parameter values of the mean field approach, the exclusion regime causes the species z 
to vanish in one time step (At = 1). This implies that, when the system is subject to 
the dynamical regime discussed in Section III Al the predator disappears. This behaviour 
disagrees with the results found by using the moment equations (see Section riIB|) . In order 



to remove this discrepancy between CML model and mean field approach, in the discrete 
time equations we use a much smaller value for the diffusion constant, namely D = 10^'^. 
By this way, we obtain a slowdown of the diffusion dynamics and, as a consequence, the 
survivance of the predator in the coexistence/exclusion dynamical regime. In order to get 
the time behaviour of the I''* and 2"''^ order moments within the scheme of the CML model, 
we consider a square lattice with = 100 x 100, using for j3{t) the time behaviour given 
in Fig. [21 Afterwards, at each time step n we calculate, from Eqs. (!26l) . (1271) . (!28|) . the new 
values of x^^^j , and the moments according to Eqs. fl33l) . (134]) . fl35l) . By iterating this 

procedure, we obtain the time series shown in Figs. El El The 1** and 2"^^ order moments 
calculated within the formalism of the CML model can be compared with the same quantities 
obtained in the mean field approach (see Figs. HI [5]). We note that the two set of time series 
are in a good qualitative agreement. According to the results obtained in the formalism of 
the moment equations, the mean values of the three species show time oscillations, whose 
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FIG. 6: In panels (a), (b) and (c) we show, respectively, the mean values, < x >'^'^\ < y >^"'\ 
< z ' , the variances, uari , vary , uar^ , and the covariances, covxy , covxy , covxy lor (7 = 0. 
The same quantities are shown in panels (d), (e) and (f) for a = 10^^^. The time series are obtained 
within the formalism of the CML model (see Eqs. (126^ . ()27p . (I28p ). The diffusion coefficient is 
D = 10^^, and 7 = 26.5. The initial values of the species concentrations are x^^^ = yf'j = 0.1, 
zf'j = 2.0 for all the sites The values of the other parameters are the same of Fig. HI A = 3, 

= 0.06, u = l,a = 0.02. 

amplitude is larger for the predator (panels a, d of Figs. El [7]). In the absence of noise, the 
2nd moments remain equal to zero (see panel b of Fig. E]), recovering the conditions 

of homogeneous distributions obtained for a = in the mean field approach (see panel b of 
Fig. H]). In the presence of multiplicative noise, no modifications occur in the time series of 
the I'** order moments (see left side panels in Figs. [6], [7]). However, for cr 7^ a symmetry 
breaking is introduced, with non- vanishing variances that are responsible for inhomogeneous 
distributions of the three species. For higher levels of the noise intensity, the amplitude of 
the oscillations remains constant in the time series of var^, vary and varz (see panels b, e 
of Figs. O [7]). These results show some difference with those obtained in the formalism of 
moment equations, where higher noise intensities cause the oscillation amplitudes of fi2oo, 
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FIG. 7: In panels (a), (b) and (c) we show, respectively, the mean values, < x >^"-', < y >^"\ 

< z >^ the variances, varx , vary , varl , and the covariances, covxy , covxy , covxy tor 
a = 10~^. The same quantities are shown in panels (d), (e) and (f) for a = 10^^. The time series 
are obtained within the formalism of the CML model (see Eqs. (126p . ()27p . (|28p ). The values of the 
other parameters and the initial conditions are the same of Fig. [6l 

/io2o and /ioo2 to become larger (see panels b, e of Figs. HI [5]). Finally, we find that for 
0" 7^ 0, temporal oscillations also appear in the time series of /ino- This agrees with the 
results of the mean field approach, revealing the presence of an anticorrelated dynamics 
between the two preys. On the other hand, /xioi and /ion remain equal to zero also in the 
presence of multiplicative noise. This behaviour, in agreement with that obtained in the 
mean field formalism, indicates that the spatial distribution of the predator is uncorrelated 
with those of each prey considered separately, but depends on the total density of preys. 
The comparison between the two approaches shows that the mean values < x >^'^\ < y >^"'\ 

< 2; >*^"^ and those obtained within the formalism of the moment equations oscillate around 
different values. Moreover, the amplitudes of the oscillations in the 2"°' order moments 
appear significantly larger in the CML model. This discrepancies can be explained recalling 
that: i) in the two approaches the stationary values are different (see Eqs. fl25l) and Eqs. (|29|) - 
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(130|) ): ii) in the mean field formalism the interaction between sites is extended to the whole 
spatial domain, conversely in the CML model the species interaction is restricted to the 
nearest neighbors; iii) the dynamics of the CML model is faster since an unitary time step 
{At = 1) is taken, instead of the time step dt = 10~^ used in the moment equations. 



IV. CONCLUSIONS 

We report a study on the stochastic dynamics of an ecosystem with three interacting 
species (two preys and one predator), described by generalized Lotka-Volterra equations. 
After considering the single site dynamics of the ecosystem, we consider a spatially extended 
domain (two-dimensional lattice) by introducing "long range" diffusive terms (diffusion oc- 
curs among each site and all the other ones). The study is performed by a mean field 
approach, in the formalism of the moment equations. The system is affected by the pres- 
ence of two noise sources, namely a multiplicative white noise and a correlated dichotomous 
noise. The role of the correlated dichotomous noise is to control the dynamical regime of 
the ecosystem (see Fig. [2]), while the multiplicative noise is responsible for the anticorrelated 
behavior of the species concentrations. The mean field approach in Gaussian approxima- 
tion enables us to obtain the time series of the 1*** and 2"'^ order moments. We compare 
the results obtained within the mean field approach with the time series calculated by a 
coupled map lattice (CML) model. The agreement is quite good, even if some discrepancies 
are present, due to the discrete nature of the CML model and the limited extension of the 
diffusive interaction (nearest neighbors) among different sites of the coupled map lattice. 
Our theoretical results could explain the time evolution of populations in real ecosystems 
whose dynamics is strictly dependent on random fluctuations, always present in natural 



environment 
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